Analysis of the Dynamics of Liquid Aluminium: Recurrent Relation Approach 
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By use of the recurrent relation approach (RRA) we study the microscopic dynamics of liquid 
aluminium at T = 973 K and develop a theoretical model which satisfies all the corresponding sum 
rules. The investigation covers the inelastic features as well as the crossover of our theory into the 
hydrodynamical and the free-particle regimes. A comparison between our theoretical results with 
those following from a generalized hydrodynamical approach is also presented. In addition to this 
we report the results of our molecular dynamics simulations for liquid aluminium, which are also 
discussed and compared to experimental data. The received results reveal that (i) the microscopical 
dynamics of density fluctuations is defined mainly by the first four even frequency moments of the 
dynamic structure factor, and (ii) the inherent relation of the high-frequency collective excitations 
observed in experimental spectra of dynamic structure factor S(k, to) with the two-, three- and 
four-particle correlations. 
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I. INTRODUCTION 

The study of dynamical processes related to atomic 
motions in liquid systems present an attractive re- 
search topic, both of theoretical and experimental origin. 
This being so, the extent of corresponding experimen- 
tal knowledge about microscopic processes in liquids has 
continuously increased and revealed new features for vari- 
ous classes of liquids, which in turn call for an appropriate 
theoretical explanation [l[ . As an intriguing example we 
mention the well-defined oscillatory modes observed in 
liquid metals occurring outside the hydrodynamic region 
and clearly detectable as corresponding high-frequency 
peaks in the experimentally relevant observable - the dy- 
namic structure factor S(k,u>) @, d, 0, IE @| • Thus, spec- 
tra of the dynamic structure factor can display, even on 
microscopical spatial scales, the triple-peak form known 
in the hydrodynamics as the 'Brillouin triplet'. 

At present various theories do exist (see, for instance, 
references @, 1, 1 E3 EH, 03, 0, El ) , which suggest an 
explanation of the 'Brillouin triplet' evolving at the finite 
values of wave-number k. Their description may be found 
in excellent reviews 0, EH EB| , and we can assess that a 
great number of these theories are based on the general- 
ization of equations of linearized hydrodynamics. How- 
ever, it is necessary to point out that hydrodynamical 
equations are a physically correct description for a large 
wavelengths and time scales, where slow processes are 
strongly pronounced and consequently provide the sig- 
nificant contribution. Therefore, any wave-number and 
especially time (frequency) extension of various thermo- 
dynamic and transport quantities or the inclusion of the 
additional non-hydrodynamical contributions should be 
invoked with some care. A need to do so is particularly 
necessary for the description of fast processes, which are a 
major component of microscopic dynamics of liquids and 



in fact may contribute appreciably to the short-time be- 
havior of corresponding relaxation functions (and/or into 
high-frequency scenario of their power spectra). The ex- 
istence of finite frequency moments of S(k,u>) can serve 
as a peculiar criterion to test and estimate the validity 
of a corresponding theoretical extension [l7| . 

In the present work we study the microscopic dynamics 
of liquid aluminium on the basis of the theory developed 
within the framework of the so-called recurrent relation 
approach (RRA) [l7|. According to the RRA, the time 
evolution of a dynamical variable depends only on basis 
vectors, which span the underlying Hilbert space of the 
studied dynamical variable and the dimensionality of the 
embedding space S [l8[ . Relaxation functions and terms 
defining the interactions appear here as projections of dy- 
namical variables on the corresponding orthogonal basis 
vectors. If the dynamical variable is a density fluctuation 
then the RRA allows one to obtain the corresponding re- 
laxation function in terms of frequency moments of the 
dynamic structure factor. 

As the specific physical system under consideration we 
choose the dynamics of liquid aluminium. This choice 
is related, foremost, to the availability of high precision 
data for the Inelastic X-ray Scattering (IXS) for liquid 
aluminium near its melting point [19J. These data are 
generally consistent with the results of molecular dynam- 
ics simulations with various models for the inter-particle 
potentials [2(| HH Nevertheless, it is necessary to 
note that the full agreement between experimental and 
simulation methods in characterization of dynamical fea- 
tures of liquid aluminium does not achieved (see p. 916 of 
Ref. (J). 

The organization of the paper is as follows: In the next 
Section we present the theoretical fundamentals; in par- 
ticular we review within our context the basic notions of 
the RRA [T3, El, El as these apply to the density fluctu- 
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ation case considered in the given work. In section HTTl the 
study of various dynamical regimes and their features for 
the case of liquid aluminium is presented. We also report 
in this section details and results of our molecular dy- 
namics simulations, which are compared with IXS-data. 
A resume with some concluding remarks is given in sec- 
tion nvi 



ii. 



THEORETICAL FRAMEWORK 



Let us consider a system of N identical classical par- 
ticles of mass m evolving in the volume V. We take the 
density fluctuations 



A(°)( k ) = J= y\ 



(1) 



as an 'initial' dynamical variable; k is the wave-vector. 
The time evolution of A"(k) is defined by the Heisen- 
berg equation 



<L4<°>(k,i) 



dt 



= i[H,A^(k,t)] = iCA^(k,t), (2) 



where C is the Liouville operator, which is Hermitian, 
TC is the Hamiltonian of the system and [H, . . .Ms the 
Poisson bracket. According to the RRA scheme jl7l . Il8j 
the quantity A^ (k, i) can be considered as a vector in 
an a priori chosen d-dimensional space S, i.e. 

d-l 

A (a \k,t) = ^2a^{k,t)f v (k), k = | k| is fixed. (3) 

u=0 

Here the 'basis vectors' fo(k), fi(k), fd-i(k) spanning 
iS are orthogonal, i.e., 



(U(k),U(k)) = (U(k),U(k))S^ 



(4) 



and are interrelated by the following recurrent relation 
(RR-I): 

U +1 {k) = itf v {k) + ^u{k)f v -i{k), ^>o, 



AAk) " (/„-!(*), /„-!(*))' 

/_i(fe) = 0, A (fc) = l. 



(5) 



Here a,u (k,t) is the time dependent projection of 
A(°'(k,t) on the v-th basis vector f v (k)\ the brackets 
(...,...) denotes the scalar product of Kubo in the embed- 
ding space S [ItJ • Then, the relaxation function between 
two variables X and Y is defined by 



(X,Y) = - (exp(XH)YUxp(-XH)X)d\, (6) 
P Jo 

where X, Y C S, (3 = (k B T)- x , k B and T are the Boltz- 
mann constant and temperature, respectively [? ]. The 



angular brackets denote the average over the canonical 
ensemble with temperature T — Note that in 

the classical case (/? — > 0, H — > 0), the relaxation function 
proportional to the usual correlation function [l7| 



(X,Y)= (XY*). 



(7) 



On the basis of the set of vectors {f n (k)} one can con- 
struct a set of dynamical variables {A^ n ^(fc)} which we 
need for the description of the evolution of the system, 
obeying 

A^(k,t = 0) = f n (k), n = 0, 1, d-l. (8) 

The corresponding time evolution by analogy with equa- 
tion ((3|) can be defined as 



d-i 



A^(k,t) = J2a^(k,t)f u (k). 



(9) 



Note that the set of functions a„(k,t) (i.e. for n = v) 
in equations §5§ and ^ are defined by projecting the dy- 
namical variable A^ n > (k, t) onto the corresponding basis 
f n {k), i.e. 



aW(M) = 



(A^(k,t),f n (k)) 
(/„(*),/«(*)) 

(#>fef),AWfeO)) 
(A(™)(fc,0),v4(")(fc,0))' 



(10) 



where the last equality is obtained by taking into account 
equation ([8]). So, one can easily see that these functions 
are normalized, i.e. 



4")(M = o) = i. 



(ii) 



Here and afterwards we simplify notation whenever the 
quantity involves the same upper and lower index by set- 
ting: d n \k,t) = a { n\k,t). 

The functions a„ (k, t) for v > n, being defined by the 
projection of the vector A^ (k, t) on the basis f v {k) with 
n 7^ v, are related to interactions between A^(k) and 
A^"\k), Further, the set a'l! l \k,t) obeys the following 
properties: 



if v > n then 
if v < n then 



a^\k,t = 0) = 0, 
aW(fc,f) = 0. 



(12) 



It is useful to note that the functions at, • (fc, t) are 
interrelated by the second recurrent relation (RR-II) [ijj , 
reading: 



A v+1 (k)a^i(k,t) 



da ( u\k,t) (n) 
J t + (13) 

n = 0, 1, 2, ... d-l, 
v = n, n + 1, ... d — 1, 
i/ > n. 
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Thus, the set of recurrent relations Q and (j 13[) yield the 
time evolution of the dynamical variable A^ n '(k), which 
occurs in the space spanned by the orthogonal basis vec- 
tors f n (k), fn+l{k), ... . 

Similar to the well-known Zwanzig-Mori formalism, 
these recurrent relations obey the equation [23[: 



j t A (n) (k,t) = A^ +1 \k,t) - A n+1 {k) 



(14) 



x / a {n+1 \k,t')A {n \k,t-t')dt', 



which constitutes the exact reformulation of equation 
(PI). The last equation can be rewritten in terms of the 
Laplace transform f(z) = e~ zt f(t)dt of the functions 

a< nJ (t) as: 



a {n \k,z) = [z + A ?1+1 (fc)a (n+1) (fc,z)]- 1 . 



(15) 



Then, the Laplace transform of relaxation function of 
the density fluctuation a^°'(k,z) can be rep resented in 
the form of a continued fraction, i.e., 24, 25, 26] 



a (Q) (k,z) 



1 



z + 



Ai(fc) 



(16) 



A a (fc) 



A 3 (fc) 



Thus, the time behavior of density relaxation func- 
tion a,(°'(k,t) is fully defined by the dimensionality d of 
the space S and by the corresponding set of parameters 
A„+i(fc) (n = 0, 1, ...). To determine A (0) (fc,i) [and 
w-°'{k,t)\ it is necessary to know the whole set of the 
parameters A„ + i(fc). 

The parameters A n+ i(fc) can be expressed in terms of 
normalized frequency moments of the dynamical struc- 
ture factor S(k,w) @,[13|, i.e., 



w y ' dtp 



t=o 



J x oo S(k,uj)dw 



(17) 



which is related, in turn, with a^ oS> (k,t) through the 
Fourier transform as 



S{k,u) 



S(k) 
2tt 



e ZUJt a^{k,t)dt. 



(18) 



Because S(k,ui) is even function of u for a classical sys- 
tem, all odd frequency moments are equal to zero. More- 
over, from equations (|13p and (|17[) and taking into ac- 



,(io) 



count equations (fTT]) and (fT2|) one finds the relations 

J 2 \k) = Ai(fc), (19) 
w (4) (fc) = A?(fc) + A 1 (fc)A 2 (fc), 
cu^(k) = A 1 (k)[A 1 (k) + A 2 (k)] 2 + A 1 (k)A 2 (k)A 3 (k), 
J 8 \k) = A 1 (fc){[A 1 (fc) + A 2 (fc)] 3 + 2A 2 (fc)A 3 (fc) 
x [Ai(fc) + Aa(fc)] + A 2 (fc)A 2 (fc)} 
+A 1 (fc)A 2 (fc)A 3 (fc)A 4 (fc), 
A 1 (fc){Ai(fe)[Ai(fe) + A 2 (fc)] 3 
+A 1 (fc)A 2 (fc)[A 1 (fc) + A 2 (fc)] 
x[Ax(fc) + A 2 (fc) + A 3 (fc)] 
+A 2 (fc)A 3 (fc)[A 1 (fc) + A 2 (fc) + A 3 (fc)] 
x [Ai(fc) + A 2 (fc) + Aa(fc) + A 4 (fc)] 
+A 2 (fc)[A 1 (fc) + A 2 (fc) + A 3 (fc)] 2 
+A 1 (fc)A 2 (fc)A 3 (fc)[A 1 (fc) + A 2 (fc)] 
+A 2 (fc)A 3 (fc)A 4 (fc)[A 1 (fc) + A 2 (fc) 

+A 3 (fc) + A 4 (fc) + A 5 (fc)]}, 



On the other hand, in the classical case, the frequency 
parameters A n (k) can be obtained in accordance with 
the help of equations (JSJl and ([7]) to read: 



(l/il 2 ) 



(l/il 2 ) 



fc R T fc 2 



(l/ol 2 



m S"(fc) ' 
fc 2 , (|/ | 2 )=5(fc), 



(20) 



By analogy with equation (|2T))) it follows that 



A a (fc) = 



k B T, 



3- 



5(fc) 



(21) 



— / V 2 u(r)[l - exp(zk • r)]g(r)d 3 r, 

?71 J 



A 3 (fc) 




A 2 (fc) 



[Ai(fc)+A 2 (fc)r 



(22) 



Here, p denotes the number density, S{k) is the static 
structure factor, g(r) is the pair distribution function, 
u{r) is the inter-particle potential and the suffix I denotes 
the component parallel to k, whereas the term !F(k) de- 
notes the combination of integral expressions containing 
the inter-particle potential with two- and three-particle 
distribution functions. In the general case, the param- 
eters A„(fc) at large nth order also contain the distri- 
bution functions of n, (n — 1), ... and n = 2. As a 
result, we can see that if the studied system is character- 
ized by strongly pronounced potential interactions then 
the problem of finding the time (frequency) dependence 
of the dynamical variables A^ n \k) and/or the functions 
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a^ n \k,t) \a^ n \k,z), see equations (|T5|) and (fT6|) ] is re- 
duced to the problem of truncating the chain of coupled 
n-th particle distribution functions [2Sj ] . 



III. DYNAMICAL REGIMES AND THEIR 
FEATURES 

A. Short wavelength dynamics 

Let us consider the spatio-temporal regime of a mono- 
atomic liquid, for which one can neglect the interaction 
between particles. Obviously, this case corresponds to 
the short-time dynamics which is restricted to length 
scales smaller than the mean free path (i.e., the regime 
of high k- values). Then, the terms characterizing the 
strength of interaction between particles are negligible, 
and the static structure factor S(k) — > 1. Upon observ- 
ing that the second (integral) term in equation (|2"Tj) and 
J-(k) in equation (f2"2")l are negligible one finds 



A 2 (fc) 



k B T, 2 



2Ai(fc), A 3 (fc) = 3Ai(fe). 



(23) 



From the first equality of equation (fTTf and equation (fl9l) 
one derives that equations (|23p represent the short-time 
behavior of a^°\k,t) described by the Gaussian function, 
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(24) 



Then, from equation (|13p one finds that other functions 

al (k, t) defining interactions of density fluctuations with 
other dynamical variables assume the following form 

af(M) = V AlW /2 ' » = 0, 1, ... ■ (25) 

Note in this case that the values of time scales deter- 
mined by ~ l/A n+ i(fc) with the increase of n become 
smaller, and the ratio between neighboring terms con- 
verges from 2 towards 1 upon increasing n; explicitly we 
have 



A„ +2 (fc) = n + 2 
A n+1 (k) n+1" 



(26) 



As a result, for this regime of large fc-values (opposite 
to the hydrodynamic limit) one can conclude that al- 
though the time scales of relaxation functions upon in- 
creasing of n also increase, they eventually become equal, 



i.e lim r 



A„+2 (k)/A n+1 (k) = 1. 



B. The fluid dynamics at intermediate wavelengths 

Obviously, the simplifications used in the previous case 
are incorrect for microscopic spatial ranges on the scale of 
order of some inter-particle distances, where interaction 



Figure 1: Main plots: the scattering intensity I(k, w) of liq- 
uid aluminium at temperature T = 973 K for different wave- 
numbers. The solid lines present the molecular dynamics re- 
sults convoluted with the experimental resolution and involv- 
ing the detailed balance condition, open circles with error 
bars are the IXS-data of reference [lSj] • Insets: the difference 
between molecular dynamics results and interpolated experi- 
mental data at the fixed k, i.e. lMD(k, w) — Iixs(k, oj), in the 
units 10~ 3 ps. 



between particles is considerable. This implies a more 
complex numerical (theoretical) evaluation of the param- 
eters A n+ i(fc), starting with the second parameter, i.e at 
n = 1 and higher. 

On the other hand, the parameters A„(fc) can be found 
from molecular dynamics simulations. This is very con- 
venient and useful tool, but it requires a priori the de- 
tailed information about the interaction of particles in 
the study system. According to this method, one can 
find the 'initial' dynamical variable (k) from simula- 
tion data [see equation ([T])] , and then receive numerically 
the structural parameters A n (k) from RR-I [see equation 
([5])]. Obviously, for higher precision the carried out re- 
sults should be averaged over time iterations. 
We have performed such numerical analysis based on 
our molecular dynamics simulations for liquid aluminium 
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Table I: Values of the static structure factor S(k), the first unnormalized frequency moment win(fc) and the scattering intensity 
at the zeroth frequency I(k,u = 0) for liquid aluminium. IXS denotes results of reference flijT ] obtained from Inelastic X-ray 
Scattering data, whereas MD indicates outcomes of our molecular dynamics simulations. 



7, / „ ,^ . — 1 \ 

k (nm J 




b[k) 




LOun (k) 




r ( 7- . . 
1 [K, U! 


= 0) (10 ps) 


TV O 


MJJ 


1~ T>7>.2 7,™ 

fcs-Z k /m 


TV C 


1VJJJ 


TVG 
JAB 


MJJ 


4.2 


0.010 


0.029 


5.288 


5.294 


5.443 


0.250 


0.805 


5.4 


0.013 


0.032 


8.741 


8.732 


11.010 


0.320 


0.826 


7.8 


0.016 


0.034 


18.239 


18.241 


23.170 


0.350 


0.837 


9.0 


0.020 


0.037 


24.283 


24.192 


32.758 


0.390 


0.844 


10.2 


0.025 


0.037 


31.190 


30.925 


40.351 


0.425 


0.896 


11.4 


0.027 


0.045 


38.961 


39.021 


47.598 


0.434 


0.910 


12.6 


0.030 


0.050 


47.595 


47.631 


59.238 


0.520 


1.174 


13.8 


0.036 


0.059 


57.092 


57.137 


72.048 


0.601 


1.253 



with the particle density n = 0.0528 A 3 at the temper- 
ature T = 1000 K. The interaction of N = 4000 particles 
imbedded in the cubic cell (L = 42.32 A) with the peri- 
odic boundary conditions was realized using the so-called 
'glue' potential 29]. The time step At applied in the in- 
tegration of the equation of motion was 10 -14 s. After 
the bringing of the system to equilibrium state, 100000 
time steps were made [30|, HH . 

Results for the static structure factor S(k) = (|/o(fc)|) 
and the first unnormalized frequency moment ujun(k) = 
(\fi(k)\) deduced from molecular dynamics simulations 
are presented in table [TJ where these quantities are com- 
pared with the corresponding values S(k) = J S(k,uj)duj 

and ujim(k) — J u> 2 S(k,u!)duj obtained on the basis of 
experimental IXS-data for the dynamical structure fac- 
tor S(k,uj) of liquid aluminium [l9| . In this table we 
also give the exact theoretical prescription for the first 

(2) 

unnormalized frequency moment uiun(k) = kBT/m. 

(2) 

As can be seen, the theoretical values of u4 n (fc) are 
in agreement with the IXS-data of reference [l9[ that is 
a quite expected, because the normalization of experi- 
mental data is performed to the first sum rules (for de- 
tails, see also reference [l[), whereas the values of uiun (k) 
obtained from simulations are overestimated in compari- 
son to IXS-data for all wave- numbers. The values of the 
static structure factor S(k) from molecular dynamics are 
also higher of IXS-data (S MD {k)/S IXS {k) ~ 1.5 - 2.9). 
As a result, the ratio between the values of the first 
structural parameter Ai(Ar) = (\.fi(k)\)/S(k) from IXS- 
data [32| and the corresponding values from molecular 
dynamics simulations is ~ 1.2 — 2.8. The analysis per- 
formed by us reveals that the difference of high-order 
structural parameters obtained from IXS-data and from 
molecular dynamics simulations is also observed, namely, 
Af D (k)/Al xs (k) ~ 2.2-6.7 and Af D (k) / A J 3 X s {k) ~ 
1.8 — 6.1. Figure [T] compares the experimental scattering 
intensity I(k, u) of liquid aluminium near melting tem- 
perature [19j | and the results of molecular dynamics sim- 
ulations. Although results of simulations are in a good 



agreement with the IXS-data for the high-frequency re- 
gion, simulation data give higher values of central peak 
in comparison to the experimental ones. The exact val- 
ues of I(k, u) = 0) from numerical simulations and from 
IXS are presented in table Q] On this basis, one can 
conclude that although the 'glue' potential of reference 
[29l ] describes correctly some equilibrium characteristics 
of liquid aluminium, it is unsuitable to reproduce the mi- 
croscopical dynamics features of liquid aluminium, and 
it is unfit for finding of the structural parameters A„(fc). 

There exists also another scheme which is based on 
a phcnomenological estimation of the structural param- 
eters A n (k). The recent experimental IXS-data of the 
intensity I(k,u>) for liquid aluminium at T — 973 K (l9j 
allows one to estimate the experimental frequency mo- 
ments and, as a result, to identify A n+ \(k). For this 
procedure we employed as a fitting function the model 
for the classical dynamic structure factor S(k,u>) as de- 
tailed in [33|, i.e., 



S(k,u) 



S(k) 



(1 - a)r 2 



sech 



tt(uj + U}q)t 2 



sech 



tt(uj - uj )t 2 



(27) 



with the model parameters a, t\, t 2 and ljq. It is worth- 
while to note that a similar approach has been used also 
in [341 ] . where the experimental S(k,uj) of liquid metals 
was fitted by a sum of three Gaussian functions in order 
to derive the first and second order memory functions 
numerically. However, as it was established during the 
fitting procedure for liquid aluminium, the model func- 
tion of 33] allows one to realize a better adjustment to 
the experimental data (see figure 1) than the combina- 
tion of Gaussian functions. As a result, after taking into 
account the detailed balance condition 

0,0, i-e., 



S q (k,u) 



hLo/ksT 

\ _ c -huj/k B T 



(28) 
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Table II: Ratio of neighboring parameters S n (k) — 
A n +i(k)/A n {k) (n = 1, 2, 3 and 4) deduced according to 
equation (|19[) on the basis of experimental frequency moments 
for liquid aluminium at T = 973 K. 



k (rnn" 1 ) 


Si(k) 


S 2 (k) 


S s (k) 


5 4 {k) 


4.2 


2.008 


1.357 


3.787 


0.935 


5.4 


2.043 


1.484 


2.839 


0.926 


7.8 


2.105 


1.344 


2.455 


0.915 


9.0 


2.466 


1.301 


2.310 


0.934 


10.2 


2.890 


1.376 


1.926 


0.964 


11.4 


2.866 


1.504 


1.680 


0.920 


12.6 


2.581 


1.465 


1.723 


0.920 


13.8 


2.519 


1.239 


1.929 


0.957 



the convolution of S q (k 1 uj) with the resolution function 
R(k, ui) was deduced to read: 



I{k,oj)=E{k) I R(k,u-u)')S q {k,u')duj, 



(29) 



wherein E(k) denotes the form factor. The resulting 
function Ifik^uj) has been adjusted to the experimental 
scattering function so that the the first two sum rules are 
obeyed identically [for S(k) and Ai(fc)]. The subscript / 
used here for the intensity I(k, u>) denotes the fact that it 
is the result of a fitting procedure. The comparison of the 
theoretical normalized second frequency parameter from 
equation (|2"D|) with the obtained one from integration of 
the fitting dynamic structure factor S(k,cj) according to 
equation (jTTJ) is depicted in the main part of figure [3] Fi- 
nally, the remaining parameters A n+ i(k) (n = 1, 2, . . .) 
were deduced from S(k,oj), which yields the best fitting 
of If(k,u>) to the experimental data. The found values 
of the first five parameters vs. k are given in the inset of 
figure [31 

During this procedure the following features for the 
studied wave-number region were identified: 

(i) although the frequency moments and parameters 
A n +\(k) are sensitive to the form of S(k,to), the ratio 
of the neighboring parameters, S n (k) = A n+ i(k) / A n (k) 
(n = 1, 2, . . .), is invariable for the various fitted func- 



tions, which reproduce the experimental data within the 
acceptable bounds; 

(ii) with increasing n for the studied wave-number re- 
gion the following feature of these parameters is observed: 
A n +i(k) > A n (k) (5 n (k) > 1). This feature is not valid 
for n = 4; 

(iii) A 4 (fc) is slightly larger than A5(fc), 0.915 < 
S 4 (k) < 0.964 (see table HJ. 

On this basis, we further suppose that all ratios of higher 
order beginning with n = 4 are approximately equal to 
1, i.e. S n (k) ~ 1 for n > 4. This suffices to obtain 
the exact equation for dynamic structure factor. From 
the physical point of view, taking into account that the 
frequency parameters can be considered as characteris- 
tics of time scales for corresponding quantities [H, [H[ , 
such an approximation means that the average time 
scales of the relaxation processes related to the dynam- 
ical variables (n > 4) are approximately equal, i.e. 
A„ 1 ^ 2 (fc) = A^( 2 (fc) for n > 4. Therefore, the time 
scales of 'neighbor' relaxation processes in this spatial 
range become equal at low levels and become strongly 
pronounced in comparison with the short wavelength dy- 
namics. As a result, in the framework of the RRA one 
can exactly obtain the relaxation function a 3 (k,t) in 
the following form: 



a^(k,t) 



1 



j 1 (2VMfe)t), 



(30) 



where Ji is the Bessel function of the first order. Fur- 
ther, relaxation functions of other dynamical variables 
and functions characterizing their interactions can be ex- 
pressed by means of integral equations. However, the 
term of the greatest interest a priori is a^°'(k,t), the 
Fourier transform of which, S(k,u>), presents a experi- 
mentally measurable quantity. Therefore, from equation 
([311)1 we obtain 



S (3) (fc,z) 



-z+ [z 2 +4A 4 (fc)] 
2A 4 (fc) 



1/2 



(31) 



Then, from equations (| 16[) and (fT3[) we can deduce the 
dynamic structure factor to read: 



2tt 



A 2 (k) 
A 3 (k) 



S(k,u) _ A 1 (fc)A 2 (fc)A 3 (fc) 
S(k) 



[4A 4 (fc) 



.211/2 



A 4 (fc)-A 3 (fc) to e + A!(k)uj 4 + A 2 (k)uj 2 

2Ai(fc), 



(32) 



A2(fc)-A 2 (fc)[2A 4 (fc)-A 3 (fc)] 



A 4 (fc)-A 3 (fc) 
A|(fc)A 4 (fc) - 2Ai(fc)A§(fc) + A 1 (fc)A 2 (fc)[2A 4 (fc) 



A 3 (fc)] 



A 4 (fc)- A 3 (fc) 



A 2 (fc), 



A 2 (fc)Al(fc) 
A 4 (fc)- A 3 (fc)' 
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Figure 2: The scattering intensity I(k,u) of liquid aluminium at temperature T — 973 K. The solid lines depict the results 
of the theoretical model (|32|l . whereas the open circles present the IXS-data of reference [jg| . The theoretical line-shapes 
have been modified to account for the quantum mechanical detailed balance condition and have been broadened for the finite 
experimental resolution effects as described in the text. 



Note that the denominator of the well known hydro- 
dynamic dynamic structure factor S(k,u>) is also repre- 
sented as the bicubic polynomial (in the variable u>) of 
equation (|32p (see reference {36}), whereas the numera- 
tor of the hydrodynamic model contains the biquadratic 
polynomial. The point is that at the crossover into the 
hydrodynamical region (as well as into the region of free- 
particle motion) the equation for the dynamic structure 
factor ([52)) becomes modified. This is so because the rules 
of ratio for the frequency parameters are also changed. 
As will be presented in the next section, the so developed 
theory yields the correct hydrodynamical asymptotic be- 
havior. 

The comparison of the theoretical results of the dy- 
namic structure factor S(k,uj) of liquid aluminium at 
T = 973 K according to equation (|32j) with data of IXS 
are presented with figure O As can be deduced from this 
figure, the theoretical results reproduce well the experi- 
mental data. 



C. High-frequency modes 

From the analysis of the experiments it is known that 
inelastic features of scattering spectra are pronounced in 
the normalized longitudinal current relaxation function 



Gj(k,t) = 



(J L (fc,0),J L (fc,i)) 



(33) 



(J L (k,o),J L (k,o)y 

which is related to the dynamic structure factor by 

S(k)Ax(k)Gj(k,oj) = uj 2 S{k,oj). (34) 
The last equation can be readily obtained from 

<9V°)(fc,i) 



A 1 (k)G J (k,t) 



dt 2 



(35) 



The quantity Gj{k,ui) possesses one minimum at lu = 
and two high-frequency maxima, being contrary to the 
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Figure 3: Main: the solid line corresponds to the theoretical results of the second frequency moment from equation (|20p . whereas 
the full circles present the values extracted by means of an integration of the resolution de-convoluted, classical S(k, ui) used 
to reproduce the experimental IXS-data (see in the text). The dashed and the dotted linespresent the hydrodynamical limit 
(see equation (|42ap with the isothermal sound velocity co(fc) = 4700 m/s, from reference [19( P and the free-particle situations 
(A/ (k) = l/vpm k), respectively. Inset: the values of the first five parameters obtained from the frequency moments of the 
dynamic structure factor S(k,u>) corresponding to the best fit of If(k, ui) to the experimental data, as detailed in the text. 



dynamic structure factor. Moreover, the side peaks in 
Gj{k,ui) correspond to those of S(k,u)). 

The inelastic features of the spectrum Gj(k, ui) are de- 
fined by the solution in regard to z = z(k) of the following 
equation (see equation (2.54) in Ref. |15l|): 



(36) 



where a^ 2 \k, z), according to the above obtained results, 
has the form 



S (2) (fc,z) 



2A 4 (fc) 



r(2A 4 (fc) - A 3 (fc)) + A 3 {k)^z 2 +4A 4 (k) 



(37) 



Generally, equation (j36|) possesses complex solutions z — 
Re[z(k))+ilm[z(k)], where Im[z(fc)] defines the positions 
of the inelastic peaks in Gj(k,ui), whereas Re[z(fc)] char- 
acterizes the widths of these peaks. 

Introducing for convenience the notation Q(k) for the 
ratio between the parameters A 4 (fc) and A3(fc), and the 
dimensionless quantity £(fc), with the latter characteriz- 
ing the frequency region: 



A 3 (fc) 



Ai(k)' 



(38a) 



(38b) 



the condition for the existence of side peaks in Eq. (|36l) 
can be rewritten in the following form (see Ref. [38l , [39 

S3): 



z 2 + zA 2 (k)- 



l + Q(k) 



= + Ai (k) = 0. (39) 

zQ(k) + v/z 2 +4A 4 (fc) 

To investigate this equation further it is convenient 
to consider the following limiting situations. First, the 
region of crossover into the hydrodynamical limit can be 
characterized by |£(fc)| <§C 1. This condition allows one to 
span the region of small frequencies (i.e. the large time 
scales) . Then, the dispersion equation takes the following 
form: 



2A\ /2 {k) 
Q(k) 



Ai(fc) 



A 2 (fc)(l + Q(fe)) 



Q(k) 



2A\ /2 {k)A l {k) 
Q(k) 



0. 



(40) 



Although the exact algebraic solution of a cubic equation 
is feasible, it is not of much practical use here because 
its inherent algebraic complexity. Let us note that equa- 
tion (|40p is similar to the dispersion equation obtained 
by Mountain (see equation (15) in reference [il[). There- 
fore, following the convergent scheme for approximating 
solutions 14jJ one finds 



Zi,2(k) = ±ic s k~rk 2 , 
z 3 (k) 



(41) 



lQ{k) ' 
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where the adiabatic sound velocity c s , the sound damping 
parameter T, and the ratio of the specific heat at constant 
pressure to the specific heat at constant volume 7 = c p /c v 
read for fc — > 0: 



2;„2 



r 

7 



V7C0, limAi(fc) = c fc 
7-lAl /2 (fc) 



7 QO) ' 
A 2 (fc)[I + Q(fc)] 
Ai(A)Q(fc) ' 



(42a) 
(42b) 
(42c) 



respectively, and Co(fc) is the isothermal sound velocity. 
The equations (|4"Tj) are the results of the hydrodynamical 
Landau-Placzek theory (43j , and the real and imaginary 
parts of the first two solutions define the positions and 
widths of the Mandelshtam-Brillouin doublet. The ap- 
proximate solutions (|4ip are valid, when the difference 
between A^fc) and As(k) is sufficiently large in compar- 
ison with differences between Aa(fc), A 2 (fc) and Ai(fc) 
(44I [45( 1 . This condition corresponds to the requirement 
of the Mountain approximation procedure (p. 208, [HI ]) 
and is related to the divergence of the frequency param- 
eters in the hydrodynamical limit. 

If one considers the large-frequency region (regardless 
of fc), i.e. the opposite extreme with z 2 / A^fc) ^ 1, then 
here equation yields 



I 

CN 

O 



o 
a 

CD 

cr 

<d 

o 

<D 
i-j 
c3 

cr 
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Figure 4: The short dashed line corresponds to Ai(fc) from 
equation (|20[) ; the open triangles present the values of 7A1 (fc) , 
which possesses the limiting behavior limfe^o 7Ai(fc) = (c s fc) 2 
and thus corresponds to the squared frequency of the Brillouin 
peak from hydrodynamical prescription (7 = c p /c v = f.4 
from references [J d3|); the filled circles denote the squared 
frequency of the side peak uj1(k) in the experimental IXS- 
spectra [ljj]; the solid line presents the theoretical predictions; 
the open squares are the values of A%(k), which are obtained 
through the second and fourth normalized frequency moments 
of the resolution de-convoluted, classical dynamic structure 
function S(k, oj) extracted from the IXS-spectra [see equation 
(I19p ]. The squared frequency of the 'solid-like' excitations 
Ai(fc) + A2(fc) = w£(fc) is presented by the stars. 



which corresponds to a typical 'instantaneous' solid-like 
response [rll42T|. 

As can be seen from the dispersion equation (|3"9"|) both, 
the widths and positions of the side peaks depend on the 
first four parameters A„(fc) (n = 1, 2, 3 and 4). In figure 
|4] the square of various frequencies are compared to the 
square of the frequency of collective excitations, w^(fc), 
extracted from the experimental data The results 

of figure |4] justify that the values of the squared frequen- 
cies of the side peaks are arranged into the intermedi- 
ate region between 'hydrodynamical prescriptions' with 
7Ai(fc) and values w£(fc) corresponding to the 'solid- like' 
behavior [equation (|4"3")) ], while the results of our model 
yield the true squared frequencies of side peak, w^(fe). 
The fact, that inelastic peaks are fully defined by the first 
four frequency parameters, indicates directly that the ex- 
istence of side peak features in the studied region depend 
on both, two- and three-particle as well as four-particle 
interactions, which are taken into account by these pa- 
rameters. Although the high-frequency motions can be 
considered as a remnant of solid dynamics, the large val- 
ues of lifetimes of the underlying excitations (in com- 
parison with the solid-state ones) are distinctive for the 
dynamics of a liquid, where two-, three- and four-particle 
correlations are strongly pronounced at the length scales 
of the order of some Angstroms. 



IV. RESUME AND CONCLUDING REMARKS 

As mentioned in the Introduction, the treatment of the 
experimental data of Inelastic Neutron Scattering (INS) 
and IXS experiments can be conveniently performed in 
the framework of generalized hydrodynamics methods. 
Recall that for the hydrodynamical regime the following 
equation holds 16] 

A 2 (fc)a< 2 >(fc,i) = ^!^) + ( 7 -I)A 1 (fc)c-^ 2t , 
pm 

(44) 

where t/l is the longitudinal viscosity and Dt is the ther- 
mal diffusion. The last equation is obtained from a rea- 
soning that is related to the hydrodynamical approach 
and is correct only for the dynamical description for large 
temporal and spatial scales, i.e t — > 00 and fc — > 0. 
However, this form of relaxation function breaks down 
when addressing the short-time features. From the point 
of view of the RRA the relaxation function correspond- 
ing to the Hermitian system can not exactly assume the 
time dependence of the form (14"4"|) , because the normal- 
ized condition (fTTj) is violated. One of the simplest gen- 
eralization of the hydrodynamical approach consists in 
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the assumption about the extended (not instantaneous) 
character of the decaying viscosity term involving vari- 
ous elementary contributions. So, in references jl], [l9| 
it was shown that for the correct reproduction of the ex- 
perimental S(k, to) (including the case of fluid aluminium 
[l9() the three exponential decay terms must be included 
at least in (k, t), 

a,W(k,t)= GjWe-WO, (45) 

j—a.^i.th 

where T a {k) and T fl (k) are the relaxation times associated 
with viscosity processes, whereas T t h{k) = 1/Z>rfc 2 is the 
relaxation time scale of thermal contribution [the second 
term in equation (|44]) ]. Q a ,n,th are the weights of the cor- 
responding contributions. Although this ansatz also has 
the violation of the normalized condition (fTTj) . it can be 
used as a sufficiently good approximation to study THz 
frequency ranges (see in reference (46l|). 

Indeed, as can be seen from figure [3J the frequency 
parameter A^(k) assumes the most significant values (in 
comparison to others) and thereby it defines the shortest 
time scales 2n / ^ A 4 (k) ~ 10 _14 s. Then the condition 
|£(fc)| <C 1 (see equation I38b[) allows one to span the fre- 
quency (time) range oj < 10 14 s _1 (t > 10~ 14 s), which 
corresponds to the microscopic dynamics and is mainly 
available to IXS and INS. As a result, we find from equa- 
tion §7$ that [13 

j = l, 2, 3, 5, ... , (46) 

where Qj{k) and Tj(fc) are expressed in terms of As(fc) 
and A4(fc) (see for details reference [3]). The given equa- 
tion at j = 3 yields the same result as equation (|45]l . 
whereas the coarse approximation at j = 1 corresponds 
to the viscoelastic model. Thus, we thereby justify that 
the expansion of the function a^ 2 \k,t) into elementary 
exponential contributions can be used as a good approx- 
imation in the description of the dynamics (but only) 
at long and intermediate time scales @ bounded by the 
picosecond regime. 

With this present study we have explored the dynam- 
ics of liquid aluminium at T — 973 K on the basis of 
the RRA, when the relaxation functions of the density 
fluctuations can be exactly extracted from the known 
relations between the frequency parameters 6 n (k) — 



A„ + i(fc)/A„(fc), (n — 1, 2, . . .). In the formulation of 
the RRA the parameters A„(fc)'s constitute the struc- 
tural characteristics of the embedding Hilbert space S of 
the investigated dynamical variables. As a result, the pe- 
culiarities of the dynamics depend directly on the dimen- 
sion d of the space S and the related parameters, A„(fc)'s. 
If, in particular, S n (k) — > 0, then the realized embed- 
ding space S can be considered as the (n+ l)dimensional 
one, and the exact solution for the relaxation function 
can be found [T3, [H, [Hj]. In the opposite situation 
S n {k) — * oo, which can be observed in the hydrodynam- 
ical limit, the corresponding exact solution causes prob- 
lems due to the divergence of the frequency parameters 
A„(fc) and A„+i(fc). 

In this work we extracted phenomenologically these 
frequency parameters on the basis of the first five fre- 
quency moments of the experimental dynamic structure 
factor. The values of these five parameters increase with 
growing index n. A similar behavior occurs for large val- 
ues of k (the free-particle case) where these parameters 
approach each other at large n. In the studied wave- 
number region this approximate identity is observed at 
ri = 4 already. The assumption that this holds as well at 
higher order n allows one to develop a theoretical model, 
which turns out to be in quantitative agreement with the 
experimentally observed IXS-data for liquid aluminium 
and which satisfies all the corresponding sum rules. We 
found that the first four even frequency moments ap- 
peared in the continued fraction of RRA (see equation 
(|16p ) as well as in the Zwanzig-Mori's formalism are nec- 
essary to restore the genuine microscopic dynamics of 
liquid aluminium. Moreover, our analysis reveals that 
the high-frequency features of microscopic collective dy- 
namics in liquid aluminium depend on two-, three- and 
four-particle correlations. The given finding can be use- 
ful in detailed investigations of equilibrium features of 
fluid aluminium, at the estimation of the many-particle 
distribution functions, and in the studying of cluster phe- 
nomena of liquid aluminium [48j near the melting. 
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